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A 
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and 2000-27990 (2000) filed February 4, 2000, the contents 
of which are incorporated hereinto by reference. 

5 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 

10 

The present invention relates to a communication 
performance calculation method and a communication 
performance calculation apparatus in a mobile 
communication system, and a blocking probability 

15 calculation method and a blocking probability calculation 
apparatus in a communication system, and a recording medium. 
More particularly, the present invention relates to a 
method and apparatus for calculating performance such as 
the blocking probability (the probability of rejection of 

20 a new call to the system) in a base station, or the 
probability of degradation in the communication 
performance due to traffic concentration, and a recording 
medium associated therewith. The present invention also 
relates to a method and apparatus for calculating the 

25 blocking probability in a mobile communication system, and 
a recording medium associated therewith. 
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DESCRIPTION OF THE RELATED ART 

In a mobile communication system like a widespread 
mobile telephone system today, the entire service area is 
5 usually divided into rather small areas called cells to 
provide services . Such a system comprises as shown in Fig. 
1, for example, a plurality of base stations 102 for 
covering divided cells 106, and mobile stations 104 for 
conducting communication with the base stations 102 by 

10 establishing radio channels with the base stations. 

The radio waves transmitted from the base stations 
102 or mobile stations 104 at certain transmission power 
travel through the space with some attenuation and arrive 
at a receiving point. The attenuation the radio waves 

15 undergo usually increases with the distance between the 
transmitting site and the receiving site. Besides, the 
propagation loss varies greatly due to the surrounding 
geography and conditions of objects because the radio waves 
undergo blocking, reflection and diffraction by large 

20 buildings, mountains or hills. On the other hand, the 
receiving side requires received power higher than a 
certain level to receive and demodulate signals at a 
prescribed level of quality. Accordingly, it is very 
important for area design in mobile communication systems 

25 to cover the service area efficiently by utilizing limited 
transmission power. 

To achieve such design, a method is often used which 
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estimates radio wave propagation conditions in the service 
areas by simulating the radio wave propagation from the 
specifications of the base stations and mobile stations 
and geographic data using a computer* Such a method is 
5 described, for example, in "Cell Design System in Mobile 
Communication" by Fujii, Asakura and Yamazaki, NTT DoCoMo 
Technical Journal Vol. 1, No. 4, pp. 28-34, 1995-01, or 
in "Total Support System for Station Establishment Design" 
by Ohmatuzawa and Yamashita, NTT DoCoMo Technical Journal 

10 Vol. 4, No. 1, pp. 28-31, 1996-04. These methods divide 
the cells into smaller subdivisions, store altitude data, 
geographic data and communication traffic data of 
individual subdivisions, and calculate the signal-to- 
interference ratio (SIR) at each receiving site or traffic 

15 for each base station. These methods employ as a multiple 
access scheme, frequency division multiple access (FDMA) , 
or time division multiple access (TDMA). 

On the other hand, as for code division multiple access 
(CDMA) , not only the propagation conditions that vary the 

20 performance of the FDMA and TDMA, but also the 

communication traffic and its temporal variations have 
great effect on the performance. Japanese Patent 
Application laid-open No. 8-191481 (1996), "Call 
Admission Control Method and Apparatus", discloses a 

25 method of deciding the admission of a new call on the basis 
of an estimate of interference at a base station, in which 
it is emphasized that the interference is an important 
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factor of deciding the performance of uplink channels in 
a CDMA system. 

Furthermore, International Publication No. 
WO98/30057 "Call Admission Control Method and Mobile 
5 Station in CDMA Mobile Communication System" discloses a 
method for a mobile station to make a call admission 
decision by transmitting information on the uplink 
interference and downlink transmission power from the base 
station to mobile stations, in which it is emphasized that 

10 the total transmission power of the base station is an 
important factor of deciding the performance for the 
downlink. Such a scheme that carries out the area design 
with considering the communication traffic and its time 
variations is new to the CDMA system. 

15 The foregoing conventional schemes have a great 

problem of being inapplicable to a CDMA system without 
change because interference caused by communications 
conducted by neighboring base stations is not counted as 
the traffic. 

20 In addition, although the interference power from the 

same radio channel reused in distant sites greatly degrades 
the performance in the FDMA or TDMA system, since it is 
not counted as the traffic, it presents a problem of 
hindering accurate performance calculation. 

25 Furthermore, since the conventional scheme does not 

consider the total downlink transmission power that serves 
as an important index in the CDMA system, it has a great 
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problem of being inapplicable to the CDMA system without 
change . 

Moreover, although the total transmission power has 
a great effect on the performance in the FDMA or TDMA system 
5 with a configuration of amplifying signals in common which 
are transmitted using multiple channels, there is no 
calculation method applicable to such a configuration, 
which presents a problem of being unable to calculate its 
perf ormance . 

10 In a communication system in which many users share 

a limited number of communication channels, there arises 
an occasion on which no communication channel is assignable 
to a user. 

In a general communication system such as a fixed 
15 telephone system or mobile telephone system, many users 
share communication resources. For example, consider an 
office telephone system shared by ten employees. The 
probability that the ten employees conduct telephone 
conversations at the same time is very small, nearly zero. 
20 Thus, the number of circuits required in the office can 
be less than ten, and the employees share the limited number 
of channels, use them when necessary, and release them 
after their conversations so that other employees can use 
them. However, it sometimes takes place that no channel 
25 is available because all the channels are busy. In such 
a case, many of the present communication systems reject 
a new call, resulting in a call loss. It is preferable 
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that the number of channels be as small as possible from 
an economical point of view. An excessively small number 
of channels, however, will increase the call loss and 
complaints of the employees, or hinder smooth processing 
5 of jobs. To satisfy such conflicting requirements, the 
Erlang B formula (see, Leonard Kleinrock, "Queueing 
Systems Volume I: Theory", John Wiley & Sons, pp. 105-106, 
1975, for example) is used to implement a sufficiently 
small blocking probability of about 1% to 3%, for example, 

10 in designing the number of channels. 

Such consideration about the office telephone system 
is also applicable to the fixed telephone system and the 
mobile communication system, as well. In particular, in 
the mobile communication system, communications between 

15 the base stations and mobile stations are established using 
radio transmission, and the resources the communications 
use are radio channels . Since the frequency band 
available for the mobile communication system is generally 
limited, the resource sharing becomes a more serious issue 

20 than in the fixed telephone network that uses wired 

circuits for information transmission. As radio channel 
access schemes the mobile communication usually utilizes, 
there are frequency division multiple access (FDMA), time 
division multiple access (TDMA) or code division multiple 

25 access (CDMA) . 

In the FDMA or TDMA system, the number of radio 
channels to be assigned to the stations can be designed 
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using the Erlang B formula as in the conventional system 
because the radio frequencies available are assigned to 
the base stations in advance. In the CDMA system, however, 
since the base stations share the same radio frequency band, 
5 the conventional method is inapplicable. 

The International Publication No. WO98/30057 n Call 
Admission Control Method and Mobile Station in CDMA Mobile 
Communication System" discloses a method of making a call 
admission decision on the basis of the uplink interference 

10 power observed by the base station and transmission power 
of the base station. It, however, only describes the 
method of making a decision as to whether the call is 
acceptable or not , and cannot obtain the blocking 
probability from the traffic applied to the system. 

15 On the other hand, there is an example that formulates 

the relationship between the applied traffic and the 
blocking probability in the CDMA system. For example, a 
paper by A. M. Viterbi and A. J. Viterbi, "Erlang capacity 
of a power controlled CDMA system", IEEE J. Select. Areas 

20 Commun., Vol.11, pp. 892-900, Aug. 1993, discloses a method 
of calculating the mean and variance of the interference 
power observed by the base station from applied traffic, 
and simply calculating the blocking probability on the 
assumption that the interference power is normally 

25 distributed. The paper calculates the blocking 
probability Packing b Y the following formula. 
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P blocking ~ Q 

where E(Z') is the mean of the normalized interference, 
Var(Z') is the variance of the normalized interference, 
5 both of which are expressed as a function of the applied 
traffic. On the other hand, A is a threshold value of the 
normalized interference, and Q(x) is defined by the 
following equation. 



A-E (Z*) 



VVa r (Z') 



(1) 



Q ( X ) =f _L,e- tV2 d t (2) 

10 

Expression (1) corresponds to calculating the 
probability that the normalized interference, which is 
assumed to have the normal distribution, exceeds the 
threshold value A. 

15 In practice, however, the interference power reduces 

when the call loss is present, and accurate calculation 
of the blocking probability is impossible without 
considering the reduction in the interference power due 
to the call loss. 

20 Fig. 14 is a block diagram illustrating the 

calculation of the blocking probability by the 
conventional technique ( disclosed by the foregoing paper ) . 
The method described in the paper calculates the blocking 
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probability without considering the reduction in the 
interference power due to the call loss. Thus, it has a 
serious problem of being unable to calculate the blocking 
probability accurately. 

5 

SUMMARY OF THE INVENTION 

The present invention is implemented in view of the 
foregoing problems. Therefore, an object of the present 

10 invention is to provide a method and the like of accurately 
and simply calculating the performance in a base station 
by calculating the traffic distribution in the base station 
in a mobile communication system. 

Another object of the present invention is to provide 

15 a method and the like of simply and accurately calculating 
the blocking probability in a mobile communication system. 
The present invention modifies the load (such as applied 
traffic (applied traffic intensity), interference power, 
transmission power of the base station) in the mobile 

20 communication system by using a dummy generated within the 
system. This makes it possible to model the reduction in 
the load due to call loss which is not considered 
conventionally, and to better simulate an actual model by 
the simplified calculation of the blocking probability, 

25 thereby improving the accuracy of the blocking probability. 
In addition, the present invention determines the amount 
of change of the dummy on the basis of the blocking 
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probability calculated and the dummy generated within the 
system, and calculates the blocking probability with 
varying the dummy. This makes it possible to bring the 
blocking probability and load reduction into balance, 
5 thereby improving the accuracy of the blocking 
probability . 

To accomplish the objects of the present invention, 
in a first aspect of the present invention, there is 
provided a communication performance calculation method 

10 in a mobile communication system which includes a plurality 
of base stations and a plurality of mobile stations for 
carrying out communication with the base stations, wherein 
an area where the mobile stations are distributed is 
divided into a plurality of subdivisions , the 

15 communication performance calculation method comprising: 
a transmission power data storing step of storing 
transmission power data of the base stations corresponding 
to the subdivisions and/or of the mobile stations visiting 
the subdivisions; 

20 a traffic intensity data storing step of storing 

traffic intensity data of the subdivisions; 

a traffic calculating step of calculating a mean and 
variance of applied traffic at the base stations; and 
a communication performance calculating step of 

25 calculating communication performance from the mean and 
variance . 

Here, the traffic calculating step may comprise: 
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a first calculating step of calculating, from the 
transmission power data of the mobile stations stored in 
the transmission power data storing step, received power 
at the base stations of signals sent from the mobile 
5 stations to the base stations; and 

a second calculating step of calculating, from the 
traffic intensity data stored in the traffic intensity data 
storing step and the received power, the mean and variance 
of the applied traffic at the base stations. 
10 The traffic calculating step may comprise: 

a third calculating step of calculating the mean and 
variance of the applied traffic at the base stations from 
the transmission power data of the base stations stored 
in the transmission power data storing step, and from the 
15 traffic intensity data stored in the traffic intensity data 
storing step. 

The communication performance calculating step may 
comprise : 

a probability calculating step of calculating 
20 probability distribution from the mean and variance of the 
applied traffic; and 

a probability decision step of calculating a 
probability that the applied traffic exceeds a 
predetermined threshold value. 
25 The probability decision step may comprise a step of 

setting acceptable interference power to the base stations 
or its constant multiple as the threshold value. 
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The probability decision step may comprise a step of 
setting a sum of acceptable interference power to the base 
stations or its constant multiple and thermal noise power 
of receivers in base stations as the threshold value. 
5 The probability decision step may comprise : 

a threshold value calculating step of carrying out 
calculation using a ratio of a sum of acceptable 
interference power to the base stations or its constant 
multiple and thermal noise power of receivers in the base 
10 stations to thermal noise power of the receivers; and 

a step of setting a calculation result in the threshold 
value calculating step as the threshold value. 

The probability decision step may comprise a step of 
setting a total sum of maximum transmission powers of the 
15 base stations or its constant multiple as the threshold 
value . 

In a second aspect of the present invention, there 
is provided a computer readable recording medium storing 
a program causing a computer to execute a communication 

20 performance calculation method in a mobile communication 
system which includes a plurality of base stations and a 
plurality of mobile stations for carrying out 
communication with the base stations, wherein an area where 
the mobile stations are distributed is divided into a 

25 plurality of subdivisions, the communication performance 
calculation method comprising: 

a transmission power data storing step of storing 
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transmission power data of the base stations corresponding 
to the subdivisions and/or of the mobile stations visiting 
the subdivisions; 

a traffic intensity data storing step of storing 
5 traffic intensity data of the subdivisions; 

a traffic calculating step of calculating a mean and 
variance of applied traffic at the base stations; and 

a communication performance calculating step of 
calculating communication performance from the mean and 
10 variance . 

In a third aspect of the present invention, there is 
provided a communication performance calculation 
apparatus in a mobile communication system which includes 
a plurality of base stations and a plurality of mobile 
15 stations for carrying out communication with the base 
stations, wherein an area where the mobile stations are 
distributed is divided into a plurality of subdivisions, 
the communication performance calculation apparatus 
comprising: 

20 transmission power data storing means for storing 

transmission power data of the base stations corresponding 
to the subdivisions and/or of the mobile stations visiting 
the subdivisions; 

traffic intensity data storing means for storing 
25 traffic intensity data of the subdivisions; 

traffic calculating means for calculating a mean and 
variance of applied traffic at the base stations; and 
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communication performance calculating means for 
calculating communication performance from the mean and 
variance . 

Here, The traffic calculating means may comprise: 
5 first calculating means for calculating, from the 

transmission power data of the mobile stations stored by 
the transmission power data storing means, received power 
at the base stations of signals sent from the mobile 
stations to the base stations; and 
10 second calculating means for calculating, from the 

traffic intensity data stored by the traffic intensity data 
storing means and the received power, the mean and variance 
of the applied traffic at the base stations. 

The traffic calculating means may comprise: 
15 third calculating means for calculating the mean and 

variance of the applied traffic at the base stations from 
the transmission power data of the base stations stored 
by the transmission power data storing means, and from the 
traffic intensity data stored by the traffic intensity data 
20 storing means . 

The communication performance calculating means may 
comprise : 

probability calculating means for calculating 
probability distribution from the mean and variance of the 
25 applied traffic; and 

probability decision means for calculating a 
probability that the applied traffic exceeds a 
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predetermined threshold value. 

The probability decision means may comprise means for 
setting acceptable interference power to the base stations 
or its constant multiple as the threshold value. 
5 The probability decision means may comprise means for 

setting a sum of acceptable interference power to the base 
stations or its constant multiple and thermal noise power 
of receivers in base stations as the threshold value. 
The probability decision means may comprise: 

10 threshold value calculating means for carrying out 

calculation using a ratio of a sum of acceptable 
interference power to the base stations or its constant 
multiple and thermal noise power of receivers in the base 
stations to thermal noise power of the receivers; and 

15 means for setting a calculation result by the 

threshold value calculating means as the threshold value. 

The probability decision means may comprise means for 
setting a total sum of maximum transmission powers of the 
base stations or its constant multiple as the threshold 

20 value . 

In a fourth aspect of the present invention, there 
is provided a blocking probability calculation method in 
a mobile communication system, the blocking probability 
calculation method comprising: 
25 a load input step of inputting a load b on the mobile 

communication system ; 

a dummy generating step of generating a dummy pd; 
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a modified load calculating step of calculating a 
modified load b' from the load b using the dummy pd; and 

a blocking probability calculating step of 
calculating a blocking probability p from the modified load 
5 b' , 

wherein the dummy generating step compares the dummy 
pd with the blocking probability p # and generates a new 
dummy pd from a compared result . 

The dummy generating step may compare the dummy pd 

10 with the blocking probability p, may make the new dummy 
pd greater than the current dummy pd when p > pd, and may 
make the new dummy pd smaller than the current dummy pd 
when p < pd. 

The blocking probability calculating step may 

15 calculate from the modified load b' a plurality of blocking 
probabilities pi, p2 , . . . , pN, where N is a natural number, 
the blocking probability calculation method may further 
comprise an estimated value calculating step of 
calculating an estimated value C from the plurality of 

20 blocking probabilities pi, p2 , pN, and the dummy 

generating step may compare the dummy pd with the estimated 
value C, and may generate a new dummy pd from a compared 
result . 

The estimated value C may equal 1 - (1 - pi) X (l-p2) 
25 X ... X ( 1 - pN) , and the dummy generating step may compare 
the dummy pd with the estimated value C, may make the new 
dummy pd greater than the current dummy pd when C > pd. 
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and may make the new dummy pd smaller than the current dummy 
pd when C < pd. 

The modified load b' may equal (1 - pd) X b. 
In a fifth aspect of the present invention, there is 
5 provided a computer readable recording medium storing a 
program causing a computer to execute a blocking 
probability calculation method in a mobile communication 
system, the blocking probability calculation method 
comprising: 

10 a load input step of inputting a load b on the mobile 

communication system; 

a dummy generating step of generating a dummy pd; 
a modified load calculating step of calculating a 
modified load b' from the load b using the dummy pd; and 
15 a blocking probability calculating step of 

calculating a blocking probability p from the modified load 
b' , 

wherein the dummy generating step compares the dummy 
pd with the blocking probability p, and generates a new 
20 dummy pd from a compared result. 

In a sixth aspect of the present invention, there is 
provided a blocking probability calculation apparatus in 
a mobile communication system, the blocking probability 
calculation apparatus comprising: 
25 load input means for inputting a load b on the mobile 

communication system ; 

dummy generating means for generating a dummy pd; 



- 17 - 




modified load calculating means for calculating a 
modified load b' from the load b using the dummy pd; and 

blocking probability calculating means for 
calculating a blocking probability p from the modified load 
5 b' , 

wherein the dummy generating means compares the dummy 
pd with the blocking probability p, and generates a new 
dummy pd from a compared result . 

Here, the dummy generating means may compare the dummy 

10 pd with the blocking probability p, may make the new dummy 
pd greater than the current dummy pd when p > pd, and may 
make the new dummy pd smaller than the current dummy pd 
when p < pd. 

The blocking probability calculating means may 

15 calculate from the modified load b' a plurality of blocking 
probabilities pi, p2, • . . , pN, where N is a natural number, 
the blocking probability calculation apparatus may 
further comprise an estimated value calculating means for 
calculating an estimated value C from the plurality of 

20 blocking probabilities pi, p2, pN, and the dummy 

generating means may compare the dummy pd with the 
estimated value C, and may generate a new dummy pd from 
a compared result . 

The estimated value C may equal 1 - (1 - pi) X (l-p2) 

25 X ... X ( 1 - pN) , and the dummy generating means may compare 
the dummy pd with the estimated value C, may make the new 
dummy pd greater than the current dummy pd when C > pd. 
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and may make the new dummy pd smaller than the current dummy 
pd when C < pd. 

The modified load b' may equal (1 - pd) X b. 

The foregoing configuration makes it possible to carry 
5 out accurate and simple calculation of the performance at 
the base stations in the mobile communication system. 

In addition, it makes it possible to carry out simple 
and accurate calculation of the blocking probability in 
the communication system. 
10 The above and other objects, effects, features and 

advantages of the present invention will become more 
apparent from the following description of embodiments 
thereof taken in conjunction with the accompanying 
drawings . 
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BRIEF DESCRIPTION OF THE DRAWINGS 



Fig. ^dTs a diagram schematically showing a 
conventional mobile communication system; 

20 Fig.^/is^a block diagram showing a configuration of 

a performance calculation apparatus of the mobile 
communication system of a first embodiment in accordance 
with the presenf invention; 

Fig. is a schematic diagram illustrating a 

25 relationship between an area divided into a plurality of 
subdivisions and records on a storage corresponding to each 
subdivision in the first embodiment; 
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Fig. 4 is a flowchart illustrating an operation when 
the present^invention is applied to an uplink channel; 

Fig. JSris a flowchart illustrating an operation when 
the present imsefition is applied to a downlink channel; 
5 Fig. 6^is a block diagram showing a configuration of 

an uplink channel performance calculation apparatus to 
which the present invention is applied; 

Fxg^f^s a block diagram showing a configuration of 
a downlink channel performance calculation apparatus to 
10 which the presprft invention is applied; 

Fig. fi^is a block diagram showing a hardware 

configuration of a blocking probability calculation 

apparatus of a second embodiment in accordance with the 

present invention ; 
Jr 

15 Fig^ST is a block diagram illustrating functions of 

the blocking probability calculation apparatus of the 



second embodijpdht ; 

Fig. \0 is a flowchart illustrating a blocking 
probability calculation processing by the blocking 
20 probability calculation apparatus of the second 
embodiment ; 

Fig. JtL is a block diagram illustrating an example 
in which a plurality of call loss factors are present; 
Fig. 12^is a block diagram illustrating functions of 
25 the blocking probability calculation apparatus of a third 
embodiment in accordance with the present invention; 
Fig. l^^is a flowchart illustrating a blocking 
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probability calculation processing by the blocking 



probability calculation apparatus of the third 
embodiment ; apci 




5 



Fig. Lfis a block diagram illustrating a conventional 
blocking probability calculation technique. 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 



The embodiments in accordance with the present 
10 invention will now be described in detail with reference 
to accompanying drawings . 

[FIRST EMBODIMENT] 



15 configuration of a performance calculation apparatus of 
a first embodiment in accordance with the present invention. 
A performance calculation apparatus 200 of the present 
embodiment comprises a main controller 202, a storage 204, 
an input/output controller 206, an input unit 208, a 

20 display unit 210 and an output unit 212, and can consist 
of a personal computer, for example. 

The left-hand side of Fig. 3 illustrates that each 
cell is divided into a plurality of subdivisions in the 
present invention. The right-hand side of Fig. 3 

25 illustrates that the storage such as a fixed disk or memory 
stores for each subdivision at least one of the 
transmission powers (Pbi and Pmi) of a base station and 



Fig. 2 is a block diagram showing a hardware 
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mobile station when the mobile station visited the 
subdivision, and traffic intensity Ai of the present 
subdivision. To calculate the received power at a 
receiving site from the stored transmission power, it is 
5 calculated using such factors stored in the storage as the 
heights above the ground of antennas of the base stations 
and mobile stations, used radio frequencies, and the 
directivities of antennas of the base stations and mobile 
stations. The received power can be calculated using an 

10 empirical formula of the received power calculation 

disclosed, for example, in M.Hata, "Empirical Formula for 
Propagation Loss in Land Mobile Radio Services " , IEEE Trans . 
on Vehic. Technol . , Vol. VT-29, No. 3, Aug. 1980. 

Fig. 4 is a flowchart illustrating an operation when 

15 the present invention is applied to an uplink channel. On 
the basis of the transmission power of the mobile station 
stored for each subdivision, the received power Ri at each 
base station, which corresponds to the signal transmitted 
from the mobile station present in the ith subdivision and 

20 received by the base station, is calculated. A step of 
calculating the mean and variance of the applied traffic 
calculates the mean of the applied traffic at the base 
station by the following equation. 

25 Average = CO + 2 {CI X Ri X Ai} (3) 

where Ai is the traffic intensity of the ith subdivision. 
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CO is a given constant, and CI is a non-negative constant. 

On the other hand, the variance of the applied traffic 
at the base station is calculated by the following 
equation . 

5 

Variance = C2 + S{C1 X ci X r± X Ri X Ai> (4) 

where C2 is a positive constant. 

Since the present embodiment is configured such that 

10 it calculates the variance besides the mean from the static 
data stored for each subdivision, it can calculate the 
degradation in communication performance and the blocking 
probability due to temporal variations in the traffic. 
More specifically, a variable k for counting the base 

15 stations is initialized to one (step S402) . Subsequently, 
a variable i for counting the subdivisions is initialized 
to one, and variables Sm and Sv for representing the mean 
and variance are initialized to zero (step S404). The 
received power of the base station is calculated for each 

20 subdivision (step S406). The calculation method is 

determined in advance as described above, in which the 
factors are used such as the heights above the ground of 
the antennas of the base station and mobile stations, the 
used radio frequencies and the directivities of antennas 

25 of the base station and mobile stations, and the foregoing 
empirical formula of the received power is used. 
Subsequently, the mean Sm and variance Sv are updated using 
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the received power (step S408). The variable i is 
incremented by one (step S410), and the same operation is 
iterated if i is equal to or less than the number of the 
intended subdivisions (step S412). If i exceeds the 
5 number of the subdivisions, the performance calculation 
is started using the calculated Sm and Sv, and a 
predetermined threshold value (step S414). The 
calculation can be carried out in various methods depending 
on performance specifications. For example, when 
10 approximating the traffic distribution by a normal 

distribution, the performance P can be calculated by the 
following equation using a complementary error function 
Erfc( 1 ). 



where T is the threshold value which can take one of the 
following values, for example: First, the interference 
power acceptable to the base station or its constant 
multiple; second, the sum of interference power acceptable 

20 to the base station or its constant multiple and the thermal 
noise power of the receiver of the base station; and third, 
a value C(l - 1/ 7] ) , where V is the ratio of the sum I to 
the thermal noise power N of the base station receiver, 
C is a given positive constant, and I is the sum of 

25 interference power acceptable to the base station or its 



P=-E r f c 
2 



[T-Sml 



(5) 



15 
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constant multiple and the thermal noise power of the base 
station receiver* Subsequently, the variable k is 
incremented by one to proceed to the calculation of the 
next base station (step S416). If k is equal to or less 
5 than the number of the intended base stations, the same 
processing is iterated for the next base station (step 
S418). If k exceeds the number of the intended base 
stations, the processing is completed. Although the 
constants CO and C2 are assumed to be zero in the present 
10 flowchart, they can be any constant. 

Fig. 5 is a flowchart illustrating an operation when 
the present invention is applied to a downlink channel. 
In this case, the mean of the applied traffic at the base 
station is calculated by the following equation. 

15 

Mean = CO + 2 {CI X Pi X Ai} (6) 

where Pi is transmission power of the base station. The 
variance of the applied traffic at the base station is 
20 calculated by the following equation. 

Variance = C2 + £{C1 X CI X Pi X Pi X Ai} (7) 

Since the present embodiment is configured such that 
25 it calculates the variance besides the mean from the static 
data stored for each subdivision, it can calculate the 
degradation in communication performance and the blocking 
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probability due to temporal variations in the traffic. 

More specifically, a variable k for counting the base 
stations is initialized to one (step S502) . Subsequently, 
a variable i for counting the subdivisions is initialized 
5 to one, and variables Sm and Sv for representing the mean 
and variance are initialized to zero (step S504). 
Subsequently, the mean Sm and variance Sv are updated using 
the transmission power stored for each subdivision (step 
S506). The variable i is incremented by one (step S508), 

10 and the same operation is iterated if i is equal to or less 
than the number of the intended subdivisions (step S510) . 
If i exceeds the number of the subdivisions, the 
performance calculation is started using the calculated 
Sm and Sv, and the predetermined threshold value (step 

15 S512). The calculation can be carried out in various 
methods depending on performance specifications. For 
example, when approximating the traffic distribution by 
a normal distribution, the performance P can be calculated 
by the following equation using the complementary error 

20 function Erfc( ' ). 



where T is the threshold value which can take the total 
maximum transmission power of the base stations or its 
25 constant multiple, for example. Subsequently, the 
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variable k is incremented by one to proceed to the 
calculation of the next base station (step S514). If k 
is equal to or less than the number of the intended base 
stations, the same processing is iterated for the next base 
station (step S516). If k exceeds the number of the base 
stations, the processing is completed. Although the 
constants CO and C2 are assumed to be zero in the present 
flowchart, they can be any constant. 

Fig. 6 is a block diagram showing a hardware 
configuration of a performance calculation apparatus 600 
of an uplink channel, to which the present invention is 
applied. In this figure, a subdivision data holding unit 
602 holds the data of the transmission power Pi and traffic 
intensity Ai of each subdivision. A propagation 
associated data holding unit 604 holds the heights above 
the ground of the antennas of the base stations and mobile 
stations, radio frequencies or directivities of the 
antennas of the base stations and mobile stations, which 
are used for calculating the received power. A received 
power calculating unit 606 reads necessary data from the 
subdivision data holding unit 602 and propagation 
associated data holding unit 604, and calculates the 
received power. A traffic calculating unit 608 calculates 
the mean Sm and variance Sv of the traffic at the base 
stations from the received power Ri for each subdivision 
calculated by the received power calculating unit 606 and 
from the traffic intensity Ai . A performance calculating 
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unit 610 calculates the performance such as the blocking 
probability and the probability of the degradation in the 
communication performance using the mean Sm and variance 
Sv calculated by the traffic calculating unit 608 and the 
5 threshold value T provided by a threshold value generator 
612, and outputs them. The calculation of the performance 
using the mean Sm and variance Sv calculated by the traffic 
calculating unit 608, and the predetermined threshold 
value T can be carried out in various methods . For example , 
10 when approximating the traffic distribution by the normal 
distribution, the performance P can be calculated by the 
following equation using the complementary error function 
Erfc( 1 ). 



P=-E r f c 
2 

15 

where T is the threshold value which can take one of the 
following values, for example: First, the interference 
power acceptable to the base station or its constant 
multiple; second, the sum of interference power acceptable 
20 to the base station or its constant multiple and the thermal 
noise power of the receiver of the base station; and third, 

a value C(l - 1/ V ) , where V is the ratio of the sum I to 
the thermal noise power N of the base station receiver, 
C is a given positive constant, and I is the sum of 
25 interference power acceptable to the base station or its 
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constant multiple and the thermal noise power of the base 
station receiver. 

Fig. 7 is a block diagram showing a hardware 
configuration of a performance calculation apparatus 700 
5 of a downlink channel, to which the present invention is 
applied. In this figure, a subdivision data holding unit 
702 holds the transmission power Pi and traffic intensity 
Ai of each subdivision. A traffic calculating unit 704 
calculates the mean Sm and variance Sv of the traffic at 

10 the base stations from the transmission power Pi and 

traffic intensity Ai. A performance calculating unit 706 
calculates the performance such as the blocking 
probability and the probability of the degradation in the 
communication performance using the mean Sm and variance 

15 Sv calculated by the traffic calculating unit 704 and the 
threshold value T provided by a threshold value generator 
708, and outputs them. The calculation of the performance 
using the calculated mean Sm and variance Sv, and the 
predetermined threshold value T can be carried out in 

20 various methods. For example, when approximating the 
traffic distribution by the normal distribution , the 
performance P can be calculated by the following equation 
using the complementary error function Erfc( 1 ). 
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where T is the threshold value which can take the total 
maximum transmission power of the base stations or its 
constant multiple, for example. 

As described above, the method according to the 
present invention comprises a transmission power data 
storing step of storing transmission power data of the base 
stations corresponding to the subdivisions where the 
mobile stations are distributed and/or of the mobile 
stations visiting the subdivisions; a traffic intensity- 
data storing step of storing traffic intensity data of the 
subdivisions; a traffic calculating step of calculating 
a mean and variance of applied traffic to the base stations ; 
and a communication performance calculating step of 
calculating communication performance from the mean and 
variance. This makes it possible to carry out accurate 
and simple calculation of the communication performance. 

In addition, according to the present invention, the 
traffic calculating step comprises a first calculating 
step of calculating, from the transmission power data of 
the mobile stations stored in the transmission power data 
storing step, received power at the base stations of 
signals sent from the mobile stations to the base stations; 
and a second calculating step of calculating, from the 
traffic intensity data stored in the traffic intensity data 
storing step and the received power, the mean and variance 
of the applied traffic at the base stations, thereby not 
only counting the interference traffic caused by 
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communications carried out by the base stations of other 
subdivisions, but also counting as the traffic the 
interference power from the same radio channel reused by- 
other subdivisions . This makes it possible to carry out 
5 accurate and simple calculation of the communication 
performance . 

Furthermore, according to the present invention, the 
traffic calculating step comprises a third calculating 
step of calculating the mean and variance of the applied 

10 traffic at the base stations from the transmission power 
data of the base stations stored in the transmission power 
data storing step, and from the traffic intensity data 
stored in the traffic intensity data storing step, thereby 
taking account of the total downlink transmission power. 

15 Using the calculation method corresponding to this makes 
it possible to carry out the highly accurate and simple 
calculation of the communication performance at the base 
stations . 

In addition, according to the present invention the 
20 communication performance calculating step comprises a 
probability calculating step of calculating probability 
distribution from the mean and variance of the applied 
traffic; and a probability decision step of calculating 
a probability that the applied traffic exceeds a 
25 predetermined threshold value. This makes it possible to 
carry out accurate and simple calculation of the 
communication performance . 
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Moreover, according to the present invention, the 
probability decision step comprises a step of setting 
acceptable interference power to the base stations or its 
constant multiple as the threshold value. This makes it 
5 possible to carry out highly accurate calculation of the 
communication performance. 

In addition, according to the present invention, the 
probability decision step comprises a step of setting a 
sum of acceptable interference power to the base stations 

10 or its constant multiple and thermal noise power of 

receivers in base stations as the threshold value. This 
makes it possible to carry out highly accurate calculation 
of the communication performance. 

Furthermore, according to the present invention, the 

15 probability decision step comprises a threshold value 
calculating step of carrying out calculation using a ratio 
of a sum of acceptable interference power to the base 
stations or its constant multiple and thermal noise power 
of receivers in the base stations to thermal noise power 

20 of the receivers ; and a step of setting a calculation result 
in the threshold value calculating step as the threshold 
value . This makes it possible to carry out highly accurate 
calculation of the communication performance. 

Moreover, according to the present invention, the 

25 probability decision step comprises a step of setting a 
total sum of maximum transmission powers of the base 
stations or its constant multiple as the threshold value. 
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This makes it possible to carry out highly accurate 
calculation of the communication performance. 

[SECOND EMBODIMENT] 
5 Fig. 8 is a block diagram showing a hardware 

configuration of a blocking probability calculation 
apparatus of a second embodiment in accordance with the 
present invention. The blocking probability calculation 
apparatus of the present embodiment comprises a main 
10 controller 801, a storage 802, an input/output controller 
803, an input unit 804, a display unit 805 and an output 
unit 806, which can consist of a personal computer, for 
example . 

The main controller 801 consists of a CPU, for example, 
15 and carries out comprehensive control of the entire system. 
The main controller 801 is connected with the storage 802. 
The storage 802 can consist of a hard disk, a flexible disk, 
an optical disk or the like. The main controller 801 is 
connected, through the input /output controller 803, with 
20 the input unit 804 such as a keyboard and a mouse, with 
the display unit 805 for displaying input data, calculation 
results and the like, and with the output unit 806 such 
as a printer for outputting the calculation results. 

The main controller 801 comprises an internal memory 
25 for storing its control program such as OS (Operating 
System) , application programs such as calculating the 
blocking probability and data required, and carries out 
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the calculation of the blocking probability using these 
programs . 

Fig. 9 is a block diagram illustrating functions of 
the blocking probability calculation apparatus of the 
5 present embodiment; and Fig. 10 is a flowchart illustrating 
a blocking probability calculation processing by the 
blocking probability calculation apparatus of the present 
embodiment . 

The blocking probability calculation apparatus as 
10 shown in Fig. 9 comprises a dummy generator 901 , an inverter 

902, a multiplier 903, a blocking probability calculating 

unit 904 and a comparator 905. 

First, the blocking probability calculation 

apparatus receives the input with a load b on the 
15 communication system through the input unit 804 (step 

S1001). In the present embodiment, it receives the 

applied traffic intensity as the load on the communication 

system. 

Subsequently, the dummy generator 901 generates the 
20 dummy pd (step S1002). A modified applied traffic 
intensity b' is calculated from the applied traffic 
intensity b using the generated dummy pd (step S1003). 
Specifically, the inverter 902 converts the generated 
dummy pd to (1 - pd) , and the multiplier 903 multiplies 
25 it by the applied traffic intensity b, thereby obtaining 
the modified applied traffic intensity b' = (1 - pd) X 
b. 
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Subsequently, the blocking probability calculating 
unit 904 calculates the blocking probability p from the 
modified applied traffic intensity b r (step S1004). In 
the present embodiment , the blocking probability p is 
5 calculated using the method (formulae (1) and (2)) 
disclosed in the paper described above. 

Then, the comparator 905 compares the dummy pd with 
the blocking probability p, and terminates the processing 
when the difference between the two is within a 

10 predetermined range, that is, when the two are considered 
to be equal (step S1005), in which case, it is considered 
that the sufficiently accurate blocking probability p is 
obtained. In other words, when the dummy pd is considered 
equal to the calculated blocking probability p, the 

15 modified applied traffic intensity b' = (1 - pd) X b is 
considered equal to (1 - p) X b, which represents the 
applied traffic intensity taking account of the call loss. 
Accordingly, the blocking probability calculated from the 
modified applied traffic intensity b' can be considered 

20 accurate . The blocking probability p thus calculated is 
displayed on the display unit 805 as needed, and is output 
to the output unit 806. 

When the difference between the dummy pd and the 
calculated blocking probability p exceeds the 

25 predetermined range, the processing is returned to step 
S1002 (step S1005). At step S1002, a new dummy pd is 
generated from the compared result of the dummy pd with 
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the calculated blocking probability p . More specifically, 
a new dummy pd is made larger than the current pd if p > 
pd, and is made smaller than the current pd if p < pd. 
After that, the processing from step S1002 to step 
5 S1004 is iterated until the difference between the dummy 
pd and the blocking probability p falls within the 
predetermined range. 

Thus, the simple and accurate calculation of the 
blocking probability can be implemented in the mobile 

10 communication system . 

Although the blocking probability calculation 
processing of the present embodiment, that is, the 
processing by the dummy generator 901, inverter 902, 
multiplier 903, blocking probability calculating unit 904 

15 and comparator 905 is carried out by means of software 
(programs), it can be implemented by means of hardware. 
In addition, the program of the blocking probability 
calculation processing can be prestored in a floppy disk, 
CD-ROM or the like so that the program can be loaded onto 

20 a memory or hard disk before its execution. 

Furthermore, although the applied traffic intensity 
is input as the load on the mobile communication system 
in the present embodiment, the interference power, the 
transmission power of the base station or the like can be 

25 input instead. In addition, although the calculation 
method disclosed in the paper described above is used for 
calculating the blocking probability, other calculation 
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methods can be utilized. 

Moreover, the present invention is applicable to any 
system such as a fixed telephone communication system and 
a mobile communication system using an FDMA, TDMA or CDMA, 
5 as long as it allows a plurality of users to carry out 
communications with sharing the resources, offering 
similar advantages . 

[THIRD EMBODIMENT] 

10 The third embodiment of the blocking probability 

calculation apparatus in accordance with the present 
invention calculates the blocking probability of each 
factor when a plurality of factors are present which bring 
about a call loss, each. The hardware configuration of 

15 the present embodiment of the blocking probability 
calculation apparatus is the same as that of Fig. 8. 

Fig. 11 is a diagram showing a case where a plurality 
of call loss factors are present. When the blocking 
probabilities of switching machines 1101 and 1102 are pi 

20 and p2, respectively, the blocking probability between A 
and B across the switching machines 1101 and 1102 is 1 - 
(1 - pi) X (i - p2). Accordingly, assuming that the 
applied traffic intensity between AB is b, the modified 
applied traffic intensity b' between AB considering the 

25 call loss is expressed as (1- (1- (1 - pi) X (l-p2))) 
Xb= ((1 -pi) X (l -p2)) Xb. Incidentally, the 
modified applied traffic intensity between AC (that is. 
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the applied traffic intensity to the switching machine 
1101) and the modified applied traffic intensity between 
BC (that is, the applied traffic intensity to the switching 
machine 1102) are equal to the modified applied traffic 
5 intensity between AB. 

Thus, the control is implemented by generating the 
dummy pd, by calculating the blocking probabilities pi and 
p2 of the switching machines 1101 and 1102 from the modified 
applied traffic intensity b' = (1 - pd) X b, and by 
10 comparing an estimated value C=l - (1 - pi) X (l-p2) 
with the dummy pd to modify the dummy pd to match them. 

Fig. 12 is a block diagram illustrating the functions 
of the blocking probability calculation apparatus of the 
present embodiment, and Fig. 13 is a flowchart illustrating 
15 a blocking probability calculation processing of the 
present embodiment . 

The blocking probability calculation apparatus as 
shown in Fig. 12 comprises a dummy generator 1201, an 
inverter 1202, a multiplier 1203, blocking probability 
20 calculating units 1204-1 - 1204 -N, a comparator 1205 and 
an estimated value calculating unit 1206. 

In Fig. 13, steps S1301 - S1303 correspond to steps 
S1001 - S1003 in Fig. 10 of the second embodiment. 

At step S1304, the blocking probability calculating 
25 units 1204-1 - 1204 -N, where N is a natural number, 
calculate the blocking probabilities pl-pN from the 
modified applied traffic intensity b' = (1 - pd) X b. In 
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the present embodiment, the blocking probabilities pl- 
pN are calculated using the method (formulae (1) and (2)) 
disclosed in the paper described above. 

At step S1305, the estimated value calculating unit 
5 1206 calculates the estimated value C=l-(l-pl)X 
(1 - p2) X...X (i - pN) from the calculated blocking 
probabilities pi - pN. 

At step S1306, the comparator 905 compares the dummy 
pd with the estimated value C, and terminates the 

10 processing when the difference between the two is within 
a predetermined range, that is, when the two are considered 
to be equal, in which case, it is considered that the 
sufficiently accurate blocking probabilities pl-pN are 
obtained. The blocking probabilities pl-pN thus 

15 calculated are displayed on the display unit 805 as needed, 
and are output to the output unit 806. 

When the difference between the dummy pd and the 
estimated value C exceeds the predetermined range, the 
processing is returned to step S1302 (step S1306) . At step 

20 S1302, a new dummy pd is generated on the basis of the 
compared result of the dummy pd with the estimated value 
C. More specifically, a new dummy pd is made larger than 
the current pd if C > pd, and is made smaller than the current 
pd if C < pd. 

25 After that, the processing from step S1302 to step 

S1305 is iterated until the difference between the dummy 
pd and the estimated value C falls within the predetermined 
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range . 

Thus, the simple and accurate calculation of the 
blocking probability in the mobile communication system 
can be implemented, 
5 Although the blocking probability calculation 

processing of the present embodiment, that is, the 
processing by the dummy generator 1201, inverter 1202, 
multiplier 1203, blocking probability calculating units 
1204-1 - 1204-N, comparator 1205 and the estimated value 

10 calculating unit 1206 is carried out by means of software 
(programs), it can be implemented by means of hardware. 
In addition, the program of the blocking probability 
calculation processing can be prestored in a floppy disk, 
CD-ROM or the like so that the program can be loaded onto 

15 a memory or hard disk before its execution. 

Furthermore, although the applied traffic intensity 
is input as the load on the mobile communication system 
in the present embodiment, the interference power, the 
transmission power of the base station or the like can be 

20 input instead. In addition, although the calculation 
method disclosed in the paper described above is used for 
calculating the blocking probability, other calculation 
methods can be utilized. 

Moreover, the present invention is applicable to any 

25 system such as a fixed telephone communication system and 
a mobile communication system using an FDMA, TDMA or CDMA, 
as long as it allows a plurality of users to carry out 



- 40 - 



communications with sharing the resources, offering 
similar advantages . 

The present invention has been described in detail 
with respect to preferred embodiments, and it will now be 
5 apparent from the foregoing to those skilled in the art 
that changes and modifications may be made without 
departing from the invention in its broader aspects, and 
it is the intention, therefore, in the appended claims to 
cover all such changes and modifications as fall within 
10 the true spirit of the invention. 
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